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ABSTRACT 

Context. Circumstellar disks are expected to form early in the process that leads to the formation of a young star, during the 
collapse of the dense molecular cloud core. It is currently not well understood at what stage of the collapse the disk is formed 
or how it subsequently evolves. 

Aims. We aim to identify whether an embedded Keplerian protoplanetary disk resides in the L1489 IRS system. Given the 
amount of envelope material still present, such a disk would respresent a very young example of a protoplanetary disk. 
Methods. Using the Submillimeter Array (SMA) we have observed the HCO"*" J = 3-2 line with a resolution of about 1". At 
this resolution a protoplanetary disk with a radius of a few hundred AUs should be detectable, if present. Radiative transfer 
tools are used to model the emission from both continuum and line data. 

Results. We find that these data are consistent with theoretical models of a collapsing envelope and Keplerian circumstellar 
disk. Models reproducing both the SED and the interferometric continuum observations reveal that the disk is inclined by 40° 
which is significantly different to the surrounding envelope (74°). 

Conclusions. This misalignment of the angular momentum axes may be caused by a gradient within the angular momentum 
in the parental cloud or if L1489 IRS is a binary system rather than just a single star. In the latter case, future observations 
looking for variability at sub-arcsecond scales may be able to constrain these dynamical variations directly. However, if stars 
form from turbulent cores, the accreting material will not have a constant angular momentum axis (although the average is well 
defined and conserved) in which case it is more likely to have a misalignment of the angular momentum axes of the disk and 
the envelope. 
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1. Introduction 

Circumstellar disks constitute an integral part in the for- 
mation of low-mass protostars, being a direct result of the 
collapse of rotating molecular cloud cores and the likely 
birthplace of future planetary systems. When a molecu- 
lar cloud core collapses to form a low-mass star, its initial 
angular momentum causes material to pile up on scales 
of ^ 100 AU in a circumstellar disk ( Ca.ssen fc MoosmarJ 
198ll : lTerebev et al.lll98l lAdams et al.lll987lll988l ). Due 
to accretion processes as the young stellar object (YSO) 
evolves, the disk grows in size while the enve l ope d issi- 
pates (e.g.. lBasulll998t lYorke fc Bodenheimeil Il999t) . In 
the more evolved stages of young stellar objects such disks 
can be imaged directly at optical and near-infrared wave- 



lengths (e.g., iBurrows et all Il996t IPadgett et aP Il999l ) 



and their chemical and dynamical properties can be de- 
rived from comparison to submillimeter spectroscopic ob- 



servations (e.g.. lThi et al.ll200l[ ). Still, this latter method 
is often not unique, but requires a priori assumptions 
about the underlying disk structure. In the earlier em- 
bedded stages this method is even more complex, be- 
cause the disk is still deeply embedded in cloud ma- 
terial so that any signature of rotation in the disk it- 
self, for example, is smeared out by emission that orig- 
inates in the still collapsing envelope. By going to higher 
resolution using interferometric observations, as well as 
observing high density gas tracers, more reliable detec- 
tions of protoplane tary disks can be made as direct imag- 
ing is approached faodrfguez et al" 19981 Qi et al. 2004 : 
IJorgensen et al. 2005b ). In this paper we present such 
high-angular resolution submillimeter wavelength obser- 
vations from the Submillimeter Array (SMA) of the cm- 
bedded YSO L1489 IRS. We show how these observations 
place strong constraints on its dynamical structure and 
evolutionary status through detailed modeling. 
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L1489 IRS (IRAS 04016+2610) is an intriguing pro- 
tostar in the Taurus star forming region [d = 140 pc) 
classified as a Class I YSO accord ing to the classification 
scheme of iLada fc Wilking (|l984j ). It is still embedded in 
a large amount of envelope material in which a significant 
amount of both infall and rotation has been observed on 
scales ranging fro m a few tenths of an AU out to severa l 
thousands of AUs (|Hogerheiidell200ll : lBoogert et al.ll2002l ). 
This large degree of rotation makes this source an inter- 
esting case study for the evolution of angular momentum 
during the formation of low-mass stars, potentially linking 
the embedded (Class and I) and revealed (Class II and 
III) stages. 

In a recent study by iBrinch et al. ( 2007t ) (referred 
to as Paper I), a model of L1489 IRS was presented 
base d on data from a large single-dish molecular line sur- 
vey ( J0rgensen et al. 2004( l. This study was motivated by 
the intriguing peculiarities seen in L1489 IRS, such as the 
unusual large size and shape of the circumstellar material. 
The aim was to constrain the structure of its larger-scale 
infalling envelope and to place it in the right context of 
the canonical picture of lo w-mass star form ation (sec for 
example several reviews in lReipurth et all 2007). 

The model presented in Paper I describes a flattened 
envelope with an inspiraling velocity field, parameterized 
by the stellar mass and the (constant) angle of the field 
lines with respect to the azimuthal direction. A spherical 
temperature profile was adopted and the model did not 
explicitly contain a disk. The mass of the envelope is 0.09 
M0 , adopted from J0rgensen et al. ( 20021 ) . The use of such 
a "global" model is sufficient when working with single- 
dish data, where the emission is dominated by the emission 
from the large-scale envelope. With this description it was 
possible to accurately reproduce all the observed single- 
dish lines. 

What the study of Paper I could not address, is 
whether a rotationally dominated disk is present on scales 
of the order of hundreds of AUs, although the amount of 
rotation which is inferred by the single-dish observations 
certainly suggests that a disk should have formed. In ad- 
dition, two specific puzzles remain about the structure of 
L1489 IRS inferred on basis of the model when compared 
to other studies in the literature. First, the best fit was ob- 
tained with a central mass of 1.35 Mq , which is a very high 
value given that the luminosi ty of the star has been deter- 
mined to be only 3.7 L© (Kenvon et al. I [1991). Second, it 
was found that the best fit inclination of the system was 
74°. T his is in agreement with the result from iHogerheii de 
(|200lh where they showed that in order to reproduce the 
observed aspect ratio the incli nation cannot be les s than 
60°. However, in a recent studv. 'Ei sner et aP ( 2005h . mod- 
eled the spectral energy distribution (SED) of L1489 IRS 
and demonstrated that this required a significantly differ- 
ent systemic inclination of only 36°. 

In this paper, we try to address these issues through 
arc second scale interferometric observations of the dense 
gas tracer HCO+ J = 3 — 2 from the Submillimeter Array. 
These observations provide information on the gas dynam- 



ics on scales ~100 AU and reveal an embedded proto- 
planetary disk. We show how careful modeling of the full 
SED from near-infrared through millimeter wavelengths 
can place strong constraints on the geometry of such a 
disk. 

The outline of this paper is as follows: Sect. [5] describes 
the details of the observations and data reduction while 
Sect. |3] presents the data. Sect. [4] introduces our model in 
which we have now included a disk and we also show how 
this model can fit all available observations including our 
new SMA data. Finally, in Sect. [5] and [H] we discuss the 
implications of our model and summarize our results. 

2. Observations and data reduction 

Observations were carried out at the Submillimeter Array 
(SMA)Q located on Mauna Kea, Hawaii. Our data set 
consists of two separate tracks of measurements in dif- 
ferent array configurations. The first track was obtained 
on December 11, 2005. This track was done in the com- 
pact configuration resulting in a spatial resolution of ^ 
2.5" with projected baselines ranging between 10 and 62 
kA. A second track was obtained on November 28, 2006, 
in the extended configuration with the resolution of about 
1" and baselines between 18 and 200 kA. The resolution of 
the two configurations corresponds to linear sizes of 350 
AU and 140 AU respectively. The synthesized beam size 
of the combined track using uniform visibility weighting 
is 0.9"x 0.7" with a position angle of 78°. 

In both tracks the receiver was tuned to HCO+ J = 3- 
2 at 267.56 GHz. We used a correlator configuration with 
high spectral resolution across the line, providing a chan- 
nel width of 0.2 kms~^ over 0.104 GHz. The remainder of 
the 2 GHz bandwidth of the SMA correlator was used to 
measure the continuum. No other lines in this band are 
expected to contaminate the continuum. We did not en- 
counter any technical problems during observations and 
the weather conditions were excellent during both tracks. 
For the compact configuration track, T225 was 0.06 and for 
the extended track T225 was at 0.08. 

Mars was used as flux calibrator in the first track, and 
Uranus for the second track. The quasars 3c454.3, 3c273, 
and 3c279 were used for passband calibration for both 
tracks, while the complex gains were calibrated using the 
two quasars 3clll and 3c84 located within 18 degrees from 
L1489 IRS. The calibrators were measured every 20 min- 
utes throughout both tracks. Their fluxes have been de- 
termined to be 3.1 and 2.6 Jy for 3clll and 3c84 in the 
compact track, and 2.3 and 2.2 Jy respectively in the ex- 
tended track. For the gain calibration, a time smoothing 
scale of 0.7 hours was used which ensures that the large 
scale variations in the phase during the track are corrected 
for. The data do not show any significant small scale varia- 
tions (i.e., rapid fluctuations in the phases or amplitudes). 

^ The Submillimeter Array is a joint project between the 
Smithsonian Astrophysical Observatory and the Academia 
Sinica Institute of Astronomy and Astrophysics and is funded 
by the Smithsonian Institution and the Academia Sinica 
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Table 1. Summary of the SMA observations. 



Source 


L1489 IRS 


0,5(2000)" 


04:04:42.85, +26:18:56.3 


Frequency 


267.55762 GHz (1.12 mm) 


Integrated line intensity 


41.85 Jy beam~^ km s~^ 


Peak line intensity 


2.04 K 


Noise level (rms) 


0.35 Jy beam~^ 


Continuum flux 


36.0±2.3 mJy 


Noise level (rms) 


3.7 mJy beam~^ 


Synthesized beam size 


0.9" X 0.7" 


(Uniform weighting) 





° Coordinates are given at the position where the continuum 
emission peaks. 



The quasars appear as point sources even at the longest 
baselines, which means that phase decorrelation due to 
atmospheric turbulence is negligible. The signal-to-noise 
of the calibrators is >50 per integration. 

The data were reduced using the MIR software pack- 
age (Qi 200^1. Due to the excellent weather conditions, 
the data quality is very high and the data reduction proce- 
dure went smooth and unproblematic. All post-processing 
of the data, including combining the two visibility se ts 



1995li^ . 



was done using the MIRIAD package (jSault et al 
Relevant numbers are presented in Tabled) 

In this paper we also make use of continuum measure- 
ments of LI 489 I RS found in the literature ranging from 
the near-infra.red ( Kenvon et al.l 1993 Eisner et al.l 2005t 



Padgett et al.lll99' 



1997t 

to (sub)millimeter wavelengths 



Myers et al 



Park fc K envon"2002': ' Whitney et al 
19871 : 



kessler-Silacci et al.l l2005l) 



( Hogerheiide &: Sandelll 



Schiev en et al.'''l994VMotte & Andre 20011 



200(j;'Moria rt 

Hogerheiide et al.l Il997: .Ohashi et al.. ,199& .Saito et al 
200lt iLucas et al.ll2000[) . 



3. Results 

The two tracks of observations provide us with two sets 
of (w, t;)-points which, when combined and Fourier trans- 
formed, samples the image {x,y) plane well from scales 
of 1 to 30" (140 to 4200 AU). Emission on spatial scales 
greater than this is filtered out by the interferometer due 
to its finite shortest spacing. This filtering is accounted 
for when comparing models to the observations. 

The measured continuum emission is shown as an im- 
age in Fig. [H Previous attempts to measure the contin- 
uum at 1. 1 mm with the BI MA interferometer were not 
successful ( Hogerheiide! 200 ll ). The SMA however, reveals 
a complex and detailed, slightly elongated, structure. 

Two reconstructed images of the HCO+ emission are 
shown in Fig. [2] using the natural and uniform weight- 
ing schemes, optimizing the signal-to-noise and angular 
resolution, respectively. In this figure, the zero moment 
map is plotted as solid contour lines and the first moment 




Aa (arcsec) 



Fig. 1. Reconstructed image of the continuum emission 
at the highest resolution (0.9" x 0.7"). The contours are 
linearly spaced with la starting at 2a and dashed lines 
are negative contour levels. 



is shown as shaded contours. It is clear that both im- 
ages reveal an elongated, flat structure. In the uniformly 
weighted image there is a large amount of asymmetry in 
the structure, which is less prominent in the naturally 
weighted image. 

The velocity contours in the lower panel of Fig. [2] are 
seen to be closed around a point which is offset with 
some 3" from the peak of the emission. There is also 
clearly a gradient in the velocity field along the major 
axi s of the object. Th is feature was previously reported 
bv lHogerheijd3 (|200lf l using interferometric observations 
of HCO+ J = 1-0, and low S/N HCO+ J = 3-2 obser- 
vations, from the BIMA and OVRO arrays. The gradient 
in the velocity field coincides almost perfectly with the 
long axis of the structure. When compared to the BIMA 
HCO+ J = 3-2 observations, the SMA data give 5-10 
times better resolution. Furthermore, the BIMA data had 
to be self-calibrated using the HCO+ J — 1-0 image as a 
model and thus the resulting image was somewhat depen- 
dent of the structure of the lower excitation emission. The 
image we obtain from the SMA data is of considerably 
better quality. 

In Fig. [31 the zero moment emission contours have 
been superposed on the near-infrared scat tered light im 



^ Available at http:/ /cfa-www. harvard. edu~cqi/mircook. html 
^ Available at^http:/ /bima.astro. umd.edu/miriad/miriad. html 



age taken by the Hubble Space Telescope (jPadgett et al 
[1999). The fi gure shows that many of the details in the 
SMA image coincide with features seen in the scattered 
light image. 

The spectral map in Fig. d] shows single spectra at po- 
sitions offset by 2" from each other. Each position shown 
here is the spectrum contained within one synthesized 
beam. The center most spectrum coincides with the peak 
of the emission in Fig. [2l Here is seen a very broad, dou- 
ble peaked line. Moving outward from the center, the lines 
become single peaked and offset in velocity space with re- 
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Uniform weight 





-5 



Act (arcsec) 



Fig. 2. Zero and first moment plots of the HCO"'" J = 3-2 
emission toward L1489 IRS. Contours are linearly spaced 
by 2(7 with a clip level of 2a. a equals 1.1 and 2.0 Jy bm~^ 
km s""'^ for the natural and uniform weighting schemes 
respectively. Negative contours appear as dashed lines. 



spect to the systemic velocity. Perpendicular to the long 
axis of the structure, the line intensity falls off quickly. 
The lines shown in Fig. [3] are reconstructed using the nat- 
ural weighting scheme. For the remainder of this paper we 
will use the uniformly weighted maps where the resolution 
is optimal. 

Fig. [5] shows the emission in a position-velocity dia- 
gram, where the image cube has been sliced along the 
major axis to produce the intensity distribution along the 
velocity axis. The PV-diagram shows almost no emission 
in the second and fourth quadrant which is a strong indica- 
tor of rotation. The small amount of low velocity emission 
seen close to the center in the fourth quadrant may be 
accounted for by infalling gas. 



4. Analysis 

With the resolution provided by the SMA it is possible to 
probe the central parts of L1489 IRS on scales of ~ 100 AU, 



Aa (orcsec) 

Fig. 3. HCO+ emission (Icr) contoi irs superposed on th e 
near-infrared scattered light image (jPadgett et al. 
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Fig. 4. Spectral map of L1489 IRS. The spectra are evenly 
spaced with 2". The bandwidth in each panel is 20 kms""'^. 
The scale on the y-axis goes from to 2 Jybm"^. 



where the protoplanetary disk is expected to be present 
(see Section [Ij. Furthermore, the HCO"'" J = 3-2 transi- 
tion traces H2 densities of 10® cm^^ or more (jSchoier et al 
I2OO5,), which are expected in the inner envelope and disk 



The model of Paper I was not explicitly made to mimic 
a disk, but rather to describe the morphology of the im- 
ages on large scales. Using this model, which works well 
on a global scale, we can explore how well it fits the SMA 
spectra when extrapolated to scales unconstrained by the 
single-dish observations. 
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Fig. 6. This figure shows three spectra from the SMA observations with model spectra superposed. The two ofF- 
positions are chosen in the direction along the long axis of the object. The offsets are chosen to be a resolution 
element. Three models are also shown: The dashed line shows the unmodified model from Paper I, the full line is the 
model from Paper I, but inclined at 40°, and the dotted line shows a model where HCO+ is absent on scales smaller 
than 200 AU. 



(D 
> 




10 5 0-5 -10 

Offset along midplane (arcsec) 

Fig. 5. Position- Velocity (PV) diagram. This plot shows 
that the emission on the scales measured by the SMA is 
entirely dominated by rotation. The black curves show 
the Keplerian velocity, calculated from the dynamic mass 
(1.35 Mq) obtained in Paper I. The curves have been in- 
clination corrected with 74° for the full lines and 40° for 
the dashed lines. Contours start at la and increase with 
1(7 = 0.12 Jy beam"i. 

Fig. [6] compares the SMA observations to the predic- 
tions of this model (dashed line). For doing this compari- 
son, the model is imaged using the (w, u)-spacings from the 
observations, so that directly comparable spectra are ob- 
tained. While the spectra away from the center are fairly 
well reproduced in terms of line width, the center position 
is too wide, i.e., the model produces velocities that are too 
extreme compared to the measured velocities. 



To check that this emission does in fact originate from 
scales less than 200 AU, a model where HCO+ is com- 
pletely absent within a radius of 200 AU is also shown in 
Fig. El with the dotted line. In this model the wide wings 
disappear and only a very narrow line is left. The two off- 
positions, which lie outside of the radius of 200 AU, are 
not affected minimally introducing this cavity. 

Fig. El also shows a spectrum which is made from the 
Paper I model, but inclined at 40° (full line). This is a 
considerably better fit than the other two models (again 
the two off-positions are little affected). Whi le the obser- 
vation s of the larger scale structure (Paper I, iHogerheiide 
(|200l[ )) demonstrates that the inclination of the flattened, 
collapsing envelope is ~74°, the SMA observations sug- 
gest that a change in the inclination occurs on scales of 
100-200 AU, reflecting a dynamically different component 
there. On the other hand, since the model of Paper I was 
tuned to match the single-dish observation on scales of 
~1000 AU or more and did not explicitly take the disk 
into account it is not unexpected that it does not fully 
reproduce the SMA observations. 

4.1. Introducing a disk model 

To reproduce also the intcrfcrometric observations we have 
modified the model from Paper I on scales corresponding 
to the innermost envelope and disk: The improvements 
consists largely of two things, namely a different parame- 
terization of the density distribution and the explicit in- 
clusion of a disk. 

The description used in Paper I has a discontinuity 
for small values of r and 6: again, this is not a problem 
when working with the large scale emission, but becomes 
a problem when interfaced with a disk where a contin- 
uous transition from envelope to disk is preferred. The 
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new parame t erizat ion we use for the envelope is adopted 
from Ulrich ( 1976[ ). In this representation the envelope 
density is given by, 



Penv{r,0) = po 



-1.5 



Rrot 

cos 9 



COS t 



COS fo 



-1/2 



Rr 



„ n . 008200 , (1) 

^ 2 cos &o r J 

where is the solution of the parabolic motion 
of an infalling particle given by r/Rroticos6o — 
cos 0) / (cos ^0 sin 2^0 ) = 1, Rrot = 150 AU is the centrifu- 
gal radius of the envelope, and po is the density on the 
equatorial plane and at the centrifugal radius. Since the 
outer radius, total mass, aspect ratio, and peak density 
are held fixed, the new parameterization is qualitatively 
similar to the Plummer-like profile used in Paper I. The 
numerical value of Rrot is chosen so that the maximum 
deviation in density in the (r, 0)-plane only reaches 20% 
in the most dense parts (outside of the disk cavity) and 
has an average deviation of less than 5%. 

In addition to this envelope, a generic disk is also in- 
troduced in the model. The density of the disk is given 

by, 



Pdisk{r,9) = 



So {r/RoY 



2ttHo 



■ exp ■ 





r cos 6 






[ Ho \ 


1 



(2) 



where 9 is the angle from the axis of symmetry, Rq is the 
disk outer radius and Hq is the scale height of the disk. The 
model thus has four free parameters: disk radius, pressure 
scale height, mass (taken as a fraction of the total mass 
of the circumstellar material), and inclination. In order to 
maximize the mid- infrared flux, we use a flat disk, i.e., 
with no increase in scale height with radius and fix the 
scale height to 0.25. The outer radius is fixed to 200 AU, 
the distance from the center of the emission in Fig. [5] to 
the position of the closed velocity contour, thereby effec- 
tively reducing the number of free parameters to two. This 
outer disk radius is not a strongly constrained parameter 
since neither the SED nor the model spectra are influ- 
enced much by changes in this value. It is only possible to 
rule out the extreme cases, where the disk cither gets so 
small that it is no longer influencing the SED or where it 
gets so big that it holds a significant fraction of the total 
mass. We find that the best fit is provided by a disk mass 
of Mdisk = 4 X 10~^ M0. A cross section of the model is 
shown in Fig. [T] 

The velocity field in the envelope is similar to the one 
used in Paper I, where the velocity field was parameterized 
in terms of a central mass and an angle between the ve- 
locity vector and the azimuthal direction, so the the ratio 
of infall to rotation could be controlled by adjusting this 
angle. The best fit parameter values obtained in Paper I of 
1.35M0 for the central mass and flow lines inclined with 
15° with respect to the azimuthal direction, are used. No 
radial motions are allowed in the disk: only full Keplerian 
motion is present in the region of the model occupied by 
the disk. 
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Fig. 7. A sHce through the model of L1489 IRS: The top 
panels show the temperature on different scales and the 
bottom panels show the corresponding density profiles. 



In order to produce synthetic observations which can 
be directly compared to the data, we use numerical ra- 
diative transfer too ls. The 3D continuum radiati ve trans- 
fer code RADMC (|Dullemond fc DominikI [20041 ) is used 
to calculate the scattering function and the temperature 
structure in the analytical axisymmetric density distribu- 
tio n. In these ca l culati ons, the luminosity measurement 
bv iKenvon et al. 1 (|l993h of 3.7 Lq is used for the central 
source. The solutio n is then "ray-traced" using RADICAL 
(jDullemond fc Tur olla 2000) to produce the spectral en- 
ergy distribution and the continuum maps at the observed 
frequencies. In both RADMC and RADICAL we assume 
certain properties of the dust. We use the dust opacities 
that give the best fit to extinction measurements in dense 
cores (Pontoppidan et al., in prep). 

The same density and temperature structure are af- 
terward given as input to th e excitation and line radia- 
tive transfer code RATRAN (jHogerheiide fc van der Tak 
2000) which is used to calculate the spatial and frequency 
dependent HCO+ J = 3-2 emission. The HCO+ mod- 
els are post-processed with the MIRIAD tasks uvmodel, 
invert, clean, and restore to simulate the actual SMA ob- 
servations. 



4.2. Modeling the continuum emission 

In Fig. (HI we compare our model to the continuum obser- 
vations. We calculate (m, w)-amplitudes at 1.1 mm which 
we plot on top of the observed amplitudes. The result is 
in good agreement with the data, suggesting that the dust 
emission is well-described by our parameterization down 
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100 150 
{uv)-dlstance (KA) 



Fig. 8. The averaged {u, -y)-amplitudes of the contin- 
uum at 1 mm in both compact and extended config- 
urations. The full line show our model. The compact 
configuration visibilities cover the (u, w)-distance up to 
60 k A. The black triangle marks the total fiux at 1 
mm ( Moriartv-Schieven et al. 1994 V The histogram indi- 
cates the zero-signal expectation values. 




10.0 100.0 1000.0 10000.0 



Fig. 9. The model spectral energy distribution assuming 
inclinations of 40° (full fine) and 50° and 74° (broken 
lines) plotted on top of flux measurements from the liter- 
ature (marked by crosses) and the Spitzer/IRS spectrum 
from 2-40 microns (red line). 



to scales of ^ 100 AU. This comparison is not very de- 
pendent on the inclination, since the continuum emission 
is quite compact. For this particular figure, an inclination 
of 40° is assumed. 



Recently, lEisner et al.1 (|2005l ) modeled the SEDs of a 
number of Class I objects in Taurus including L1489 IRS. 
The model used in their work is parameterized differently 
from ours, but it esse ntially describ es a similar structure. 
As pointed out above. lEisner et al.l find a best fit inclina- 
tion of 36° in contrast to the inclination of 74° found in 
Paper I. 

To test the result of Eisner et al. ( 20051 ) we calculated 
the SED using the described model with only the system 
inclination as a free parameter (Fig. [9]). The best fit is 
found for an incli nation of about 40°, in good agreement 
with the result of Eisner et al.l ( 20051 ). The models with 
inclinations of 50° and 74°, are also plotted in Fig. [5] in 
order to show how the SED depends on inclination. 

Note that the quality of our be st fit to the SED i s 
comparable to the fit presented by lEisner et al. (|2005h , 
whereas the fit becomes rapidly worse with increasing 
inclination larger than 40°, th us bringing about dis- 
agreement with the results from faogerheiide 200 if ) and 
Paper I. The three different models plotted in Fig. [HI can- 
not be distinguished for wavelength above ~60 fim, which 
corresponds to a temperature of about 40 K using Wiens 
displacement law. This temperature occurs on radial dis- 
tances of approximately 100 AU from the central object 
which means that outside of this radius, the SED is no 
longer sensitive to the inclination. We interpret this be- 
havior as a change in the angular momentum axis on disk 
scales, causing the disk to be inclined with respect to the 
envelope. 




-6-4-2 2 



Fig. 10. Three transitions of HCO+ observed by single- 
dish telescope (see Paper I for details) with our model 
superposed. The quality of the fit is comparable to the fit 
in Paper I although a disk model has now been included. 



4.3. Modeling the HCO^ emission 

In Paper I it was found that the characteristic double peak 
feature of the IICO+ J = 4-3 line could not be reproduced 
with an inclination below 70° . On the other hand including 
a disk inclined by 40° into the model of Paper I does not 
alter the fit to the single dish lines ('Fig fTU)) : the geometry 
and the velocity field of the material at scales smaller than 
~300 AU does not influence the shape of these lines. This 
fit is not perfect though as it still overestimates the red- 
shifted wing in the J = 3-2 line and the width of the J = 
1-0 line slightly, but the quality of the flt is similar to that 
presented in Paper I. 

We need to test how well this tilted disk model works 
with the line observations from the SMA. It was shown in 
Sect. |4] that the off-position spectra are weakly dependent 
on changes in the inclination and that the same is true 
when using the model where the disk is tilted with respect 
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Fig. 11. The SMA spectrum towards the center of 
L1489 IRS with our model, consisting of a disk which is 
inclined with respect to the envelope, superposed. As in 
Fig.[Hl the model has been imaged with the (u, w)-spacings 
from the observations. 




50 100 150 200 

(uv) — distance (k\) 



Fig. 12. Similar to Fig. M but for the HCO+ J = 3-2 
measurements. The triangle marks the total flux obtained 
from the single dish observations. The amplitudes have 
been averaged over ±5 kms~^ relative to the systemic ve- 
locity. The histogram indicates the zero-signal expectation 
values. The insert is an enlargement of the longest base- 
lines. 



to the envelope. The central position however, is seen in 
Fig. [TT]to be very well reproduced in terms of line width 
and wing shape by this model. The model spectrum is 
slightly more asymmetric than the data, which means that 
the infall to rotation ratio is not quite correct on disk 
scales. The magnitude of the velocity field projected onto 
our line of sight is correct though, since the width of the 
line is well fit. 

We also compared the averaged (u, w)-amplitudes to 
evaluate the quality of our model similar to the pro- 
cedure for the continuum. Fitting the (u, w)-amplitudes 
in this way tests whether the model produces the right 
amount of emission at every scale. The best fit is shown 
in Fig. [121 The zero-spacing flux has been shown in this 
plot as well, marked by the black triangle. This method 
has previously been used to study the abundance varia- 
tions of given molecules in protostellar envelopes, in par- 
ticular imaging directly where significant freeze-out occurs 



in prot ostellar envelopes ()j0rgensen.2OO4Hj0rgensen et al. 
The model does a good job reproducing the ob- 



served HCO"'" brightness distribution on almost all scales 
with a constant abundance, except around 20-40 kA where 
it overestimates the amount of observed flux. This corre- 
spond to a radius of about 1000 AU, that is well outside 
of the disk. On the other hand single-dish observations 
of low-mass protostellar envelopes suggest that these are 
scales where significant freeze-out in particular of CO may 
occur at temperatures < 20 — 30 K, which because of 
the gas-phase chemical relation between CO and £[00+ 
also r eflects directly in the o bserved distribution of HCO"*" 
(e.g., IJorgensen et"al]l2004L lioOSb l. Despite this, the dif- 
ference between the model prediction and observations is 
small with this constant abundance, suggesting that the 
amount of freeze-out is small in L1489 IRS. This is in 



agreement with the single-dish studies of iJorgensen et al 
(|2002l . I2OO4L l2005bh . who generally found little depletion 
in the envelopes of L1489 IRS and other Class I sources 
in contrast to the more deeply embedded. Class 0, proto- 
stars. 

The model parameter values obtained in this section 
including values from Paper I are summarized in Table [21 
The best fit values that are derived in this paper are tuned 
by hand rather than systematically optimized by a min- 
imization method. Therefore we can only give an estimate 
on the uncertainties. However, Fig. [H shows that the SED 
fit gets rapidly worse when changing the inclination, and 
given the strong dependence of the envelope inclination 
on the single-dish lines, we estimate that both angles are 
accurate to within 10°. For the disk mass and radius, the 
uncertainties on the parameter values are less well defined 
for reasons given in Section 031 but we estimate an accu- 
racy within a factor of 2 in each of these parameters. 

5. Discussion 

It seems that depending on the physical size scales that 
we probe, the solution favors a different inclination. When 
taking all available observations into account we need to 
introduce a model where the angular momentum vector 
of the disk is misaligned with the angular momentum vec- 
tor of the (rotationally flattened) envelope. We adopt a 
two component model (illustrated in Fig. [T51) . disk and 
envelope, but in reality, the angular momentum axis may 
indeed be a continuous function of radius on the scales 
of the disk (< a few hundred AU). The origin of such a 
misalignment can actually be explained in a simple way 
by considering the initial conditions of the collapse. 
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Table 2. Model parameter values 



Envelope 


Value 


Outer radius*^ 


2000 AU 


Envelope mass^ 


0.093 Mq 


Inclination'^"'' 


74° 


HCO+ abundance" 


3.0x10"^ 


Disk 


Value 


Radius'* 


200 AU 


Disk mass* 


4x10"^ Mq 


Inclination'' 


40° 


Scale height'^ 


z=0.25R 


Central mass" 


1.35 M0 


Central luminosity'^ 


3.7 Lq 


Distance'^ 


140 pc 



" Values from Paper I 
Fitted values 

Values adopted from the literature 




500 AU 

Fig. 13. An (exaggerated) illustration of the proposed 
model where the angular momentum axis is changing with 
radius. The line of sight is illustrated by the dashed line. 



The state of a pre-stellar core before the gravitational 
collapse begins is typically described by a static sphere 



with a solid body rotation perturbation (iTerebev et al 



19841) . Due to strict angular momentum conservation, such 
a model would be perfectly aligned throughout the col- 
lapse and consequently disk formation models are often 
descr ibed numerically by axisym metric representations 
(e.g., Yorke fc Bodenheimer 1999f ). 

However, there is no a priori reason why a pre-stellar 
core should rotate as a solid body. Of course the cloud has 
an average angular momentum which on a global scale de- 
termines the axis of rotation. The cloud does not collapse 
instantly to form a disk tho ugh, but ra ther from the in- 
side and out as described by IShu ( 1977 ). and therefore a 
shell of material from deep within the cloud, which may 
very well have an average angular momentum that is dif- 
ferent from the global average, will collapse and form a 
disk before material from further out has had a change to 
accrete yet. Actually, if the parental cloud is turbulent. 



it is to be expected that the accreting material has ran- 
domly oriented angular momentum and a misalignment of 
the angular momentum on different scales is more likely 
than a perfect alignment. In this way, a system similar to 
the model proposed to describe L1489 IRS in this paper 
can be formed. 

To obtain a gradient in the angular momentum as ini- 
tial condition for the collapse, one may consider a cloud 
which is not spherically symmetric or one which does not 
collapse around its geometrical center but rather around 
an over dense clump offset from the center. The former sce- 
nario may indeed be true for the case of L1489 IRS which 
is seen to be dynamically connected to a neighboring cloud 
(which was also modeled in Paper I to explain the excess 
of cold emission seen in some of the low- J lines). Actually, 
the uniformly weighted moment map (Fig. [2]) shows sig- 
nificant asymmetry with a secondary emission peak to the 
north-east. 

It is interesting to note that the HCO^ emission from 
the SMA observations agree very well with the near- 
infrared image (Fig. [3]): the secondary peak, a few arc sec- 
onds northeast of the main continuum and HCO"*" peaks, 
nicely coincides with a bright spot in the scattered light 
image and also the shape of the cavity towards the south 
follows the contours of the IICO+ emission. The same de- 
tails are not revealed in the naturally weighted SMA im- 
age, in which shorter baselines are given more weight. We 
thus conclude that by going to the extended SMA config- 
uration, it is possible to probe structure in YSOs on the 
same scales as can be resolved by large near-infrared tele- 
scopes such as the Hubble Space Telescope. The structure 
seen in both images, however, also emphasizes that for 
a fully self-consistent description of L1489 IRS, a global 
non- axisymmetric model has to be considered. 

Another (non-exclusive) explanation of a misaligned 
disk would be that L1489 IRS formed as a triple stellar 
system and, that due to gravitational interaction, one of 
the stars was ejected. L1489 IRS would thus be a b inary 
system as suggested bv lHogerheiide fc Sandell (|200d ). The 
loss of angular momentum due to the ejection would result 
in a rearrangement of the remaining binary, which could 
"drag" the inner viscous disk a long. Such a scenario has 
been investigated numerically bv lLarwood et aT (1996). In 
the case of L1489 IRS, it would have to be a very close bi- 
nary with a separation of no more tha n a few AUs since the 
near-infrared scattered light image (jPadgett et al.lll999f) 
with a resolution of ^0.2" (30 AU), does not reveal mul- 
tiple sources. The ratio of binary separation to disk ra- 
dius is t herefore significantly l ower than the cases investi- 
gated bv lLarwood et al" ( 19961 ) and it is therefore not clear 
whether similar effects could be present in L1489 IRS. 

If indeed the system is binary, it would resolve the 
issue that we find a central mass of 1.35 MfT^ w hile the 



luminosity is estimated to be 3.7 Lq (jKenvon et al.lll993l ). 
A single young star that massive would require a much 
higher luminosity, but two stars of 0.6-0.7 M0 would fit 
nicely with the estimated luminosity since luminosity is 
not a linear function of mass for YSOs. 
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To test whether L1489 IRS indeed is a binary and to 
constrain the innermost disk geometry on AU scales, mea- 
siircmcnts of temporal variability in super high resolution 
are needed. The timescale for variations of a possible bi- 
nary is of the order of years depending on the exact sepa- 
ration of the stars (assuming scales of 1 AU) . Such obser- 
vations would, for example, be feasible with ALMA. 

6. Conclusion 

In this paper, we have presented high angular resolution 

interferometric observations of the low-mass Class I YSO 
L1489 IRS. The observations reveal a rotationally dom- 
inated, very structured central region with a radius of 
about 200-300 AU. We interpret this as a young Keplcrian 
disk which is still deeply embedded in envelope material. 
This conclusion is supported by a convincing fit with a 
disk model to the SED. 

We conclude further that the inclination of the disk is 
not aligned with the inclination of the flattened envelope 
structure, due to the possibility that L1489 IRS is a binary 
system and/or that the average angular momentum axis 
of the cloud is not aligned with the angular momentum 
axis of the dense core that originally collapsed to form the 
star(s) plus the disk. 

We find that a disk with a mass of 4 x 10~^ M©, a 
radius of 200 AU, and a pressure scale height of z=0.25R 
is consistent with both the SED and the HCO+ observa- 
tions. Only a small amount of chemical depletion of HCO"'" 
is allowed for, due to a slight over-estimate of the (u, v)- 
amplitudes at 20-40 kA by our model, in agreement with 
the results from previous single-dish studies and the na- 
ture of L1489 IRS as a Class I YSO. 

The combination of a detailed modeling of the SED 
with spatially resolved line observations, which contains 
information on the gas kinematics, appears to be a very 
efficient way of determining the properties of disks, es- 
pecially embedded disks which are not directly observable. 
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